On the microscopic theory of the exciton ring fragmentation 
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The description is presented for the dependence of the indirect exciton condensate density at the 
ring as a function of the polar angle at zero temperature with the involvement of the processes of 
formation and recombination of the excitons. In particular, starting from the quasi one-dimensional 
Gross-Pitaevskii equation with a spatially uniform generating term, we derive an exact analytical 
solution yielding the fragmentation of an exciton ring which is probably observed in the experiments. 

PACS numbers: 71.35.-y, 71.35.Lk, 03.75.Hh 

Keywords: Indirect excitons; Bose-Einstein condensation; Pattern formation 



1. Introduction. In the recent experiments of the 
Butov's group P, Q the macroscopically ordered coher- 
ent state (MOCS) of the indirect exciton system has been 
observed in the coupled quantum wells. The state repre- 
sents a periodic fragmentation of the ring formed in the 
intersection domain of the surplus densities of electrons 
in one well and photoexcited holes in another 0, 11, S]- 
The fragmentation is observed below some critical tem- 
perature (Tc ^ \K) and within a certain range of the 
laser pumping power when the density n of indirect ex- 
citons (lEs) at the ring is sufficiently large and the exci- 
tons are well defined (n ~ lO^cm^^, an average spacing 
between lEs at the ring is (iex-ex ^ rT^I"^ ~ O.l/im, 
the IE Bohr radius is as ~ 10~^//r7i, hex-ex/ as ^ 10). 
The typical parameters of the fragmented ring are as fol- 
lows [2|: average ring radius is i? '--^ Vfipm, the ring 
width is Ai? ^ 10/im, the average radius of a fragment 
is r/r « Ai?/2, and a number of fragments on the ring is 
Nfr ~ 50. The evidence for the MOCS coherence results 
from a clear interference pattern of the IE recombination 
radiation from one of the ring fragments Q- 

At present, there exist two different phenomenological 
theories of the MOCS formation. The first theory [3| is 
based on the assumption that the IE system represents 
a degenerate two-dimensional quasi-ideal Bose gas and 
the rate of aggregating an electron and a hole to form 
an indirect exciton is proportional to the factor (1 -I- iVo), 
where Nq is the number of the lEs in the zero momentum 
state. Under definite conditions this leads to the instabil- 
ity against fiuctuations of the IE ring density. However, 
the theory seems to be intrinsically inconsistent since, on 
the one hand, the dipole-dipole interaction between the 
excitons is neglected, what is allowable only in the low 
density limit, and, on the other hand, the exciton con- 
densate density should be sufficiently large so that the 
contribution from the stimulated processes would be es- 
sential. 

The second theory @ is based on that the total con- 
tribution of the repulsive dipole-dipole and attractive 
van der Waals interactions becomes attractive at the dis- 
tances smaller than several IE radii. According to this 
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theory, the ring fragmentation is due to formation of the 
islands of electron-hole liquid and is not associated with 
the Bose-Einstein condensation of the lEs. (We also men- 
tion paper Q in which the IE energy distribution is sug- 
gested to be the Boltzmann one and, along with an at- 
tractive pair potential, the repulsive three-body interac- 
tion between lEs is considered. At the same time the na- 
ture of the latter is not explained.) However, the results 
of the experiment unambiguously evidence for the re- 
pulsive character of the IE interaction. This is consistent 
with the idea that just this interaction results in a rapid 
IE thermalization (along with the nonconservation of a 
transverse momentum) and in the screening of the sur- 
face defects in the quantum wells plane. Note that any 
extension beyond the dipole approximation while con- 
sidering the interaction between the lEs is equivalent to 
the deviation from the Bose statistics, which should self- 
consistently be taken into account. 

As is seen from the pictures of the spatial distribution 
of the IE luminescence (see 3), there exist two differ- 
ent spatial periods on the fragmented ring, namely, the 
size of a fragment and the size of a dark region separating 
adjacent fragments. Therefore, it is clear that the den- 
sity as a function of the polar angle should be expressed 
via a periodic function with nonsymmetric half-periods, 
i.e., an eUiptic function. So, if p ((/>) is an elliptic func- 
tion, it satisfies equation (yP'^ — Pi (p), where Pi [p) 

is a quartic polynomial in p. To the same form one can 
reduce the one-dimensional stationary Gross-Pitaevskii 
equation without an external potential, which, as is well 
known, describes a condensate of Bose particles. Hence, 
the idea arises that the ring fragmentation can be ex- 
plained assuming that under experimental conditions an 
IE condensate exists on the ring. In addition, the for- 
mation and recombination of excitons on the ring can be 
taken into account by their interaction with some scalar 
generating, or "source-drain", field independent of coor- 
dinates. Formally, this is similar to the contribution from 
the dipole interaction of atoms with the electromagnetic 
field. If the number of the lEs on the ring does not vary 
in time, a concept of energy of the system is meaningful 
and one can find the stationary states of the condensate 
as well as the corresponding density spatial distributions. 
In this Letter, we suggest the simplest microscopic the- 
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ory of the MOCS based on the thing that the presence of 
BEC in the IE system (in spite of the BEC depletion due 
to repulsive exciton-exciton interaction) as well as the 
influence of the processes of IE formation and recombi- 
nation on the ring result in the realization of a stationary 
excited state of the condensate, which is not spatially uni- 
form. The indirect excitons are treated as genuine Bose 
particles. In addition, we suppose that the radial proflle 
of the condensate density is governed by the processes 
which have no direct relation to the ring fragmentation. 
This allows us to consider a quasi one-dimensional prob- 
lem. 

2. Polar angle dependence of the condensate 
density. The Hamiltonian for the exciton system with 
the source reads 



H = Hj+ I dfi'+(f,t)^if(r,t)+ 



(1) 



-J dfdr^/+{f,t)ii+{r,t)U{r- r')*(r , t)*(f, i). 
Here ^'(r, t) is a Bose field operator obeying the com- 



mutation relations 



*(fi,t),*+(r2,t) 



5 {ri - r2). 



0. For simplicity, we suggest that 

the repulsive pair interaction is local, i.e., J7(r) = XS (f) 
and A > 0. In fact, this means either the Born approx- 



imation for the potential U{r) with A = y U{r)dr, or 

the gas approximation na\jj ^ 1, where n is the two- 
dimensional IE density and a2D — mX/ (47rfi,^) is an exact 
s-scattering length [11]. The contribution of the source 
is given by 



Hj= dr 



'(j(i)*+(f,t) + J* (0*(r,t) 



(2) 



where J (t) — \ J\ e'"^*\ It is natural to assume that the 
time dependence of the source is the same as for the wave 
function of the ground condensate state in the lack of the 
source, a{t) — —fit. We suggest that the quantities |J| 
and are time- independent. 

Let us consider the equation of motion for ^-operator. 



assuming zero temperature 



of the system. Here, the most fraction of excitons is in the 
Bose-condensate state and thus we neglect nonconden- 
sate particles. Then ^'-operator becomes the c-number 
and obeys the equation 



2m 



Al^-l ^-f J, 



(3) 



which for J = turns into the Gross-Pitaevskii equa- 
tion. Supposing that the radial profile of the conden- 
sate density is unvaried, we consider a quasi ID prob- 
lem. Then, the genuine density "in amplitude" refers to 
the result obtained as ntrue/n ~ o^ex-ex/^R ~ 10"'^. 
For the polar frame with the origin at the ring center 
« R-'^d'^/dcj)^, dr K 27ri?Ai?d0 = Sd(t), where cj) is 



the polar angle, R is the ring radius, and R (n). The 
aim of the work is to find a qualitative behavior n{(j)) at 
various magnitudes of the parameters. 

If one measures energy in units En = It? / {2mR^) , time 
in Tfl = h/Eji, and density in , the constants in 
Eq.|[3l) are dimensionless 



-VI/'' + A 1*1'* + ./. 



(4) 



Let us seek for the solution ([4]) as * = pe*-^, where ~ 
n and / are the density and phase of the condensate, 
respectively. Then, for the real and imaginary parts of * 
the following equations are valid 

Pf + P (/;)' - P'U + V + \J\ cos = 0, (5) 

p + 2p'^/; + - I J| sin (a - /) = 0. (6) 

Multiplying both sides of Eq.|[6l) by 2p and integrating 
over 0, we obtain that the rate of variation of the particle 
number N equals 



AT^y" dr2pp = J dr2p\J\sm{a- f) . 



(7) 



In what follows, we are interested in the stationary cur- 
rentless states alone: N — 0, f'^ = 0. Then f = a + nl, I 
is an integer. So, 

p';^ + PP-xp^-\j\{-if ^0, (8) 

in which a (t) = —fit is taken into account. In the case of 
the homogeneous density (p^ = 0, p = po) and | J| = 0, 
/i has a meaning of the chemical potential of the system: 

P^Xpl. (9) 

Note that for | J| = even in the unhomogeneous case the 
quantity p can literally be treated as a chemical potential 
for the ground state of the system. In this case it can be 
introduced via *(i) ~ exp {—ipt). 

For J 7^ the uniform density is determined from the 
equation 



pp^-xpi-\j\{-\y ^Q. 



(10) 



Employing the Cardano formulae, one can find from 
the analysis that there exist three regions in which: 
(i) homogeneous solution is absent (even I, p < p^, = 
3A(|J|/(2A))'/^), (ii) home geneous solution is unique 
(even I, p = /i*, and odd I at p ^ /i*), and (iii) there 
exist 1 or 2 various homogeneous solutions (arbitrary I, 
P > p*). 

For the unhomogeneous case 7^ 0, the first integral 
of Eq. ([HI) is given by 

(P^j 



\Xp*-pp^ + 2\J\{-lfp + C. (11) 



Here C is a constant independent of (j). Substituting p 
po • X and integrating, we have 



dx 



v/a;4 - 2pxx'^ + 2Jxx + Cx 



±^/^Vo (<?!'-0o) = 0, 



(12) 
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where ^\ 



|J|(-i)' 
— I 



C. 



c 



xi,X2,X3, and X4 be the roots of equation — S/i^^;^ 
2JxX + Cx = 0. From HI) it follows 



Let 



^ y/Xl - X3y/X2 - X ^ 

\/Xi - X2y/X3 - x' 



-a(j>, 



where a = y^^-x^^x^-x^, - 
Here F [x, k) is the Jacobi elliptic integral of the first 
kind with modulus k. Solving the latter equation in x, 
we arrive at 

X2 {xz - xi) + xz {xi - X2) sm? (a^, fc) 
Xz — Xi + {xi — X2) srfi I a4>, k 1 

where sn{x, k) is the Jacobi elliptic sine with modu- 
lus k. (Note that the quantities xi, X2, xs, 2:4 should 
be chosen so that a; ^ 0, > and < /c^ ^ 1. 
Otherwise, this unhomogeneous solution does not exist.) 
For k — > 1, the period of s-n?{x,k) tends to infinity 
and the density is independent of the polar angle. For 
k — > 0, sn^{x,k) — > sin^(x). The position of the ex- 
trema of function x (4>) is determined by the quantity 
^ = (xi - X2) {x2 - X3) {xi - X3). li £, = 0, x'^ = 0. For 

^min ^ mi2l^ X (^P'^) = 



^ > 0, the points of minima 



Kik] 



(21 + 1), ^(^f^^'^) = xi,l 



X2 and maxima 

being an integer. For ^ < 0, in these formulae one should 
replace 4>f"^ ^ ^f^"". 

It is seen that, in the general case, the condensate den- 
sity n(0) = PqX^ {(p) is a periodic function (see Fig. 1). 
The solution obtained allows us to find the period of den- 
sity oscillations or, this is the same, the angular size of a 
fragment 



A(j>fr 



2K{k) 



(14) 



where K (k) is the complete elliptic integral. The number 
of fragments at the circumference of the ring Nfr = -^fy 
should be a positive integer. One can, in principle, find 
the constant C = C {R) from this condition. 

Note that the spatially unhomogeneous solutions are 
also possible in the lack of the source (| J| — 0). Then 



n{(t>) = {poxi) sn^ 




^pI Mx2,ko 



(15) 



where ^0 = f;, 



= Ma ± Vma - C'a 



3. Conclusions and discussion. To conclude, we 
have obtained a set of solutions. All of them, both ho- 
mogeneous and unhomogeneous, correspond to the mini- 
mum of the energy's functional of the system, i.e., to the 
stationary solutions (see Appendix). This simple theory 




Figure 1: The examples of the behavior of the condensate 
density n versus polar angle (j) for xi = 1.2, X2 = 0.8, X3 = 0.4, 
X4 = 1.21 (solid line) and for xi = 1.2, X2 = 0.75, X3 — 0.6, 
0:4 = 2.1 (dashed line). 



cannot uniquely predict what a solution from the set will 
be realized in fact. However, it is clear that the spatially 
homogeneous solution (if it exists at J 7^ 0) refers to the 
absolute minimum of energy and unhomogeneous ones 
correspond to the steady excited states of the conden- 
sate. 

In experiments [J, [3| the IE formation-recombination 
processes lead to the condensate density fiuctuations in 
the radial direction. As is shown in Ref. at zero 

temperature in the cigar-shaped condensate due to the 
parametric resonance there takes place one-way transfer 
of energy from the transverse (radial) modes to the lon- 
gitudinal ones. Closing this "cigar" in a ring, we come 
to the point that the spatially non-uniform distribution 
of the condensate density is nothing else than a standing 
wave corresponding to the macroscopically filling of the 
longitudinal modes. The consideration of the standing 
wave as a steady excited state of the condensate agrees 
with the experimental fact that the existing fragmenta- 
tion starts to vanish with decreasing the pumping power 
but for all that the ring radius would vary negligibly. 

Note that the fragmentation of the exciton ring can 
be explained in another way, assuming that the lEs are 
in the condensate state. Since the "ground" state of the 
exciton system is degenerate over the IE spin configu- 
rations, it represents a mixture of the condensates with 
different coupling constants Xij > for the mixture com- 
ponents (here i,j are the indexes of the mixture compo- 
nents). As is known [13, 14], if for i ^ j the condition 
Xfj > XiiXjj is fulfilled, the mixture components sepa- 
rate from each other in space. (On the observation of 
this effect in atomic condensates see, e.g., Ref. [l3|.) It 
is possible that such a separation for the ring geometry 
is observed in the experiments. 
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Recently, the observations similar to MOCS ones have 
been reported in Ref. ^6|. The microscopic derivation 
(starting with the electron-hole Hamiltonian) of Eq. ^ 
with J = is the same as in Ref. • 
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Appendix. The total average energy of the system 
has form 

m -^Jd<t> ((P^)' - Mp' + + 2 I J| (-1)' . 



(16) 



The minimum of the functional E[p] is determined by 
Eq.|[8]) and equals 

En.in[p] = ^Jdc^{2{p'^f-C), (17) 


where C is the constant from Eq. lfTT]) . In the case of 
the homogeneous density, involving Eq. (fTO|) . we arrive at 
£"111111 [po] = MPo ~ fAPo = —Cq. This corresponds to the 
absolute energy minimum provided that 

l.Jd^{p'^)' >^{C~-Co). (18) 
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